Mon. Not. R. Astron. Soc. 000. HHTol (20061 Printed 5 February 2008 (MN I*TeX style file v2.2) 



The redshift distribution of SWIFT Gamma-Ray Bursts: 
evidence for evolution. 
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ABSTRACT 

We predict the redshift distribution of long Gamma-Ray Bursts (GRBs) with Monte 
Carlo simulations. Our improved analysis constrains free parameters with three kinds 
of observation: (i) the log AT — logP diagram of BATSE bursts; (ii) the peak energy 
distribution of bright BATSE bursts; (iii) the HETE2 fraction of X-ray rich GRBs and 
X-ray flashes. The statistical analysis of the Monte Carlo simulation results allow us to 
carefully study the impact of the uncertainties in the GRB intrinsic properties on the 
redshift distribution. The comparison with SWIFT data then leads to the following 
conclusions. The Amati relation should be intrinsic, if observationally confirmed by 
SWIFT. The progenitor and/or the GRB properties have to evolve to reproduce the 
high mean redshift of SWIFT bursts. Our results favor an evolution of the efficiency 
of GRB production by massive stars, that would be ~ 6 — 7 times higher at z ~ 7 than 
at z ~ 2. We finally predict around 10 GRBs detected by SWIFT at redshift z > 6 
for a three year mission. These may be sufficient to open a new observational window 
over the high redshift Universe. 
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1 INTRODUCTION 

Gamma-ray bursts (hereafter GRBs) are powerful flashes 
of high-energy photons, that travel undisturbed from 
cosmological distances to earth, where they can be eas- 
ily detected. Therefore, independently of their physical 
origin, GRBs can be used to probe the distant Universe 
jLamb fc Reichaxd |20o3 iBlain fc Nataraianl l2000t iTotanil 
Il999l: IWiiers et alJ 11993) . In addition, there has recently 
been growing observatio nal evidence that lo ng GRBs occur 
in star-forming regions feloom et al.ll200lT) and that some 
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distance and their association with massive stars motivate 
two questions: (i) how is the long GRB rate related to 
the star formation rate (SFR) ? (ii) what is the expected 
detection rate of long GRBs at high redshift ? On the one 
hand, a better understanding of the first issue is a necessary 
step to be able to use GRBs to directly trace the star 
formation history in the universe. Observations give strong 
indications in favor of an association of long GRBs with 
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massive stars. However, the actual physical conditions (e.g. 
mass, metallicity, rotation, binarity) for a star to trigger 
a burst are not currently known. Thus the GRB and star 
formation rates may be related in a non trivial way. On 
the other hand, GRBs could be produced in the Universe 
in association with first stars and are thus exp ected up to 
very high redshifts (e.g. iBromm fc Loeb|l2005t) . Therefore, 
they are potentially important tools of investigation of the 
cosmic evolution; they are probably unique direct probe 
at z > 6.5. First, the afterglow emission may be used to 
study the ionization and metal enrichm ent histories of th e 
intervening intergalactic medium (e.g. IChen et alJ I2006T) . 
Then, they can give insight into galaxy formation and 
evolution: they allow us to detect very faint galaxies at high 
redshift, that othe rwise would have eluded detection (e.g. 
iBerger et al.ll2006T) . In turn, the detection of high redshift 
galaxies can better constrain the SFR where current results 
are most uncertain. Finally, they might be important to 
constrain the small -scale power spectru m of primordial 
density fluctuations jMesinger et al.ll2005l) . 

The standard approach to investigate such questions 
is to assume a GR B comoving rate, lumino sity function 
and spectrum (e.g. iPorciani fc Madaul 1200 it thereafter 
PM01) and to constrain the free parameters of the model 
by observed data, especially the GRB log N — log P diagram 



2 F. Daigne, E.M. Rossi and R. Mochkovitch 



JSchmidtl Fl999V We follow this approach, by the mean 
of Monte Carlo simulations. Such an approach allows a 
realistic parametrization of the intrinsic GRB properties, 
a more accurate treatment of detection criteria for several 
instruments and a study of the impact of the uncertainties 
in the GRB physics on t he predicted GRB rate . Com- 
pared to prev i ous studies llLamb fc Reichartl l200 o} PM01. 
[ Firmani etafl 12004 iGuetta et alJ 120051 : iNataraian et alJ 
2005), we add two more observational constraints that 
allow us to better determine the model parameters: the 
observed peak energy distribution of bright BATSE bursts 
and the fraction of X-ray rich GRBs (hereafter XRRs) and 
X-ray flashes (hereafter XRFs) in the HETE2 GRB sample. 
Together with the data obtained during SWIFT first year, 
this method allows us to address the two questions that 
motivated this work. 

The paper is organized as follows. In section |21 we de- 
scribe our assumptions on the GRB intrinsic properties (co- 
moving rate, luminosity and spectrum). Then, we describe 
our Monte Carlo simulations in sectional We especially de- 
tail our detection criteria for several instruments (BATSE, 
HETE2 and SWIFT) and the observations we use to con- 
straint the free parameters of the model. Our results are 
discussed in section^Jand sectionals the conclusion. 
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Figure 1. The three star forma tion rate s consi dered in this paper. 
Observed data are taken from iHopkinJ 120 



2 GRB INTRINSIC PROPERTIES 

Our Monte Carlo simulations are necessarily based on some 
assumptions regarding the GRB intrinsic properties. In some 
cases, the large uncertainties in these properties lead us to 
consider several scenarios. In this way, we can study the im- 
pact of this poor physical understanding of the GRB phe- 
nomenon on our predictions, especially on the GRB redshift 
distribution. 

2.1 Comoving rate 

Searching for high redshift galaxies has extended the 
number of rest fram e UV-bright systems over the interval 
^ z < 10 (e.g. iGiavalisco et al.1 liooll: iBouwens et alJ 
2005). These data have been used to probe the redshift 
evolution of the SFR. It seems established that the SFR 
density grows by an order of magnit ude from z = to 
z ~ 1 and levels off around z ~ 2 f see iHopkinsI (120041) and 
Fig. for a recent collection of data). The shape of the 
star formation history is more uncertain for z > 3. Recent 
investigations of photometric dropouts at z > 6 suggest a 
strong decline of more than an orde r of magnitude from 
z = 3 to z ~ 10 iBouwens et alJl20o3. 120051 : iFontana et all 
l2003l:IStanwav et alJl2004) . However, the assessment of this 
result is subject to the poorly constrained dust obscuration 
at such high redshifts. 

Since long GRBs are very likely related to the death 
of massive stars, the most simple assumption is to adopt 
a GRB rate proportional to the SFR. Following PM01, we 
take 

Kgrb = k x ftsN , (1) 
where TZghb is the GRB comoving rate and TZsn is the type 



II supernova comoving rate. We assume the same initial 
mass function (IMF) as in PM01. As this IMF is constant 
with time, the supernova rate is proportional to the SFR : 

■Rsn = 0.0122 Mq 1 x SFR , (2) 

the lifetime of massive stars being neglected. We first con- 
sidered for the SFR a simple fit of observed data up to z ~ 6 
jHopkinsll2004D . This case is thereafter called SFRi. As the 
true evolution of the SFR above z ~ 2 — 3 is still uncertain, 
we have also considered two alternative rates : SFR2 , where 
the SFR is constant above 2 ~ 2 and SFR3, where the SFR 
still increases above z ~ 2. In this last case, we arbitrarily 
cut the GRB rate at z max = 20. These three SFRs are plot- 
ted in Fig.[T] They are in f act ve ry close to those adopted by 
PM01. lLamb fc Reichartl J200(J) considered two cases, one 
being close to SFRi and one being close to SFR2 with an 
additional very intense com ponent at z > 6 — 7. T he pre- 
ferred SFR in the results bv lNataraian et alJ (1200 51) i s close 
to our SFR 3 up to z ~ 6 - 7. Finally iFirmani et all (12004ft 
adopted a rather different approach. The SFR in their study 
is a three free parameters function. It is adjusted together 
with the luminosity function under the joined constraint of 
the BATSE log N - log P diagram and the "observed" GRB 
redshift distribution provided by the lu minosity-variability 
relation jFenimore fc Ramirez-Ruizll200Cf) . 

The ratio of GRB explosions to type II supernovae, k, 
is a free parameter of the model. Here 7?.grb has to be un- 
derstood as the comoving rate of GRBs pointing towards 
us. The true GRB rate in the Universe (as well as the true 
GRB / SN ratio) is obtained by multiplying this rate by a 
correcting factor for the b eaming, which is of the order of 
((Q./4-k)- 1 ) ~ 500 - 1000 jFrail et alJl200lh . 



The redshift distribution of SWIFT GRBs: evidence for evolution. 



3 



io 56 

10 55 

in 

M 1054 

8 io 63 

10 52 
2.5 



ft 
_o 

CO 



1.5 



log N — log P only 



- | mi | rTTTTTTl] — I' I I II U| |-4- I 1 I I — I I | 



i nun 1 1 1 urn 1 1 H i m n 



mil i i I I I I nl 



10 47 10 48 1 49 1 50 1 61 10 52 
Lmin [erg/s] 



1 1.5 2 2.5 
Slope 



10 56 
10 55 

w 

W 10 64 

to 

i io 63 

10 52 



2.5 - 



a 

CO 



1.5 



1 - 



E p only 



I llllllj — I I llllll| — I I llllllj — I I llllllj — I I llllllj — I I I I I | I I I I | 



1 1 li nn 1 1 1 urn 1 1 nun n 



□in i ' i mill I I I I 



1 1.5 2 2.5 
Slope 



10 47 IO 48 1 49 1 50 1 51 10 52 
k,,„ [erg/s] 




Figure 2. Parameter space (SFR.3, Amati-like relation) : we plot the location in the par ameter space of the best models (see 
section Top-left: models that fit the B ATSE log N - log P diagram <Stern et al.ll20n(l Enii^ : Top-right: models that fit the peak 
energy distr ibution of bright BAT SE bursts iPreece et al.ll200d) : Bottom-left: models that fit the fraction of XRFs and XRRs observed 
by HETE2 iSakamoto et alJl200.'fl : Bottom-right: models that fit the three constraints. 



2.2 Luminosity function 

The GRB luminosity function is currently quite uncertain. 
We assume a power-law distribution, which does not vary 
with redshift: 



p(L) oc L s for L min < L < L max . (3) 



Notice that p(L) is the distribution of the isotropic equiva- 
lent luminosity of GRBs that would be obtained by a perfect 
(no threshold) detector on Earth. This is of course different 
from the true distribution of L for all GRBs in the Universe, 
which could be obtained from p(L) after correction for view- 
ing selection effects. With p(L) there are three more model 
parameters: the slope S and the minimum and maximum 
luminosities i m i n and L ma x- Previous studies also used such 



4 F. Daigne, E.M. Rossi and R. Mochkovitch 



a power-law luminosity function jLamb fc Reichard l200d) 
or a Schechter luminosity function fPMOl. TNefeiraianet^Ll 
2005). This latter is very close to a single power-law ex- 
tending to I/max = +00 with an additional low-luminosity 
tail below L m i n t hat does not contribu te much due to an 
exponential cutoff. iFirmani et alJ l|2004f ) also used a power- 
law distribution, but including the additional possibility 
that the typical GRB luminosity scales with (f + z) v with 
v ~ 0.8 — 1.2. We also briefly tested such evolutionary effect 
(see section l4"Tt . 



2.3 Spectrum 

The intrinsic photon spectrum is as sumed to be a br oken 
power-law with a break at energy E p ijBand et al.ll993T) and 
a low (resp. high) energy slope —a (resp. — /3). We checked 
that a m ore realistic spect rum shape, the so-called "GRB" 
function feand et alJll993h . has only a small impact on our 
results and does not affect our conclusions. The probability 
distributions of the low-energy and high-energy slopes fol- 
lows the obse rved distribution o f a and /3 in a set of long 
bright GRBs iPreece et al.ll2000f) . We have checked a poste- 
riori that the simulated distributions of the "observed" low- 
and high-energy slopes in long bright GRBs were very close 
to the intrinsic distributions. 

As the observed peak energy di stribution of these same 
bright burst is close to log-normal (IPreece et "ai]|200rJl . the 
most simple assumption is to adopt an intrinsic log-normal 
distribution for E p , with a dispersion a and a mean value 
Epfi. This is our first scenario (hereafter "log-normal E p dis- 
tribution"). We assume a — 0.3 dex (we checked that our 
results do not depend too much on this assumption) and we 
keep Epfl as a free parameter. On the other hand, there are 
some evidences that the intrinsic peak energy could be corre- 
lated with the intrinsic luminosity. Therefore, we have con- 
sidered a second scenario (hereafter "Amati-like relation"), 
where 

/ T \° A3 

E p = 380 keV = , (4) 

p ^1.6 x 10 52 erg s' 1 J w 

with a normal dispe rsion a = 0.2 dex, in agreement with 
the o bserved relation <|Yonetoku et al.p004: Ghirland^e^^J 
120051. see however lNakar fc Piranl2005l : lBand fc PreecefeOOa 
who tested this relation against BATSE data and concluded 
that selection effects were dominant). In this second sce- 
nario, no free parameter is introduced for the spectrum. 

Notice that the treatment of the spe ctrum in this work 
is imp roved compared to previous studies: lLamb fc Reichartl 
( 2000) adopted a single power-law spectrum of slope a = — 1 
which of course is a very poor app roximation for high 
redshift GRBs. Other studies (PM01, IFirmani et"ai]|2004t 
iNataraian et all 120051) adopt a standard GRB spectrum 
with fixed slopes a = — 1 and /3 — —2.25 and a constant 
break energy Eb in the GRB rest frame which is a priori fixed 
to 511 keV (E p = E b x (2+a)/(a-p) ~ 410 keV in the GRB 
rest frame), without any a posteriori check that the simu- 
lated obs^rveo^Egdistrihution is in agreement with BATSE 
data. IFirmani et alJ (I2004T) also consider a case including an 
intrinsic Amati-like relation of slope 0.5, but without any 
dispersion. 



Table 1. 


Best models: parameters. 




SFR 


log Z/ min log L max S 


log k 




Amati-like relation E p ex L 43 




1 


49.9 ±0.5 53.7 ±0.4 1.70 ± 0.08 


-5.4 ±0.3 


2 


50.0 ±0.5 53.7 ±0.5 1.68 ±0.10 


-5.5 ±0.3 


3 


50.3 ±0.7 53.5 ±0.4 1.54 ±0.18 


-6.0 ±0.2 


log-normal distribution peak energy distribution 


1 


50.2 ±0.9 53.6 ±0.8 1.62 ±0.18 


-5.6 ±0.3 


2 


50.2 ±1.1 53.6 ±0.9 1.62 ± 0.27 


-5.7 ±0.3 


3 


50.5 ± 1.3 53.7 ±0.9 1.52 ± 0.48 


-6.2 ±0.2 



SFR log£ p ,o 

log-normal distribution peak energy distribution 



1 2.74 ± 0.08 

2 2.73 ±0.08 

3 2.79 ±0.08 



3 MONTE CARLO SIMULATIONS 

In our Monte Carlo simulations, each generated GRB is 
given a redshift, a luminosity and a spectrum, according to 
the specific intrinsic distributions that have been described 
above. Then observed properties (observed peak energy, 
peak flux) can be computed and compared to the sensitivity 
of different instruments. In this study, we focus on BATSE, 
HETE2 and SWfFT. Our detection criteria are specified in 
section |3~T1 This allows us to derive the expected distribu- 
tions of observed redshifts, peak fluxes and peak energies 
for these instruments, for comparison with real data. It also 
allows to reconstruct the source frame properties (e.g. lu- 
minosity, peak energy) of detected bursts. In a second step, 
the free parameters of the model are adjusted to fit several 
observational constraints. This is explained in section 

3.1 Detection criteria 

We consider the detection by three instruments: BATSE, 
HETE2 and SWIFT. For BATSE, we apply the detection 
efficiency as a f unctio n of the peak flux, e(P), given by 
iKommers et al.l i2000t) . so that we can compare our results 
with the log N — log P diagram published by these authors. 
We also compare our s imulated bursts wi t h the log N — log P 
diagram published bv lStern et"aH fcOOCt 120021) . Their sam- 
ple goes farther towards low flux bursts and puts there- 
fore more constraints on the GRB rate at high redshift. 
Their published diagram is already corrected for detection 
efficiency. In addition we apply a threshold P50-300 kcv > 
5 ph cm -2 s _1 defining a sub-sample of "bright BATSE 
bursts", to compare our simulated GRBs with the res ults 
of the BATSE spectroscopic catalog iPreece et aljEoOpT) . In 
our simulations, this threshold typically corresponds to 5 % 
to 10 % of BATSE bursts, in agreement with observations. 
For HETE2, we adopt the same threshold of 1 ph cm~ 2 s _1 
both for the WXM (2-10 keV) and for FREGATE (30-400 
keV) (J.L. Atteia, private communication). 

Swift uses two detection modes: one based on a flux 



The redshift distribution of SWIFT GRBs: evidence for evolution. 



5 



Table 2. Best models: mean redshift. 



Rate All SWIFT Bright SWIFT 
Amati-like relation E p oc L° 43 



SFR1 3.1 1.6 1.3 

SFR2 8.0 1.9 1.5 

SFR3 10.5 3.3 2.1 



log-normal peak energy distribution 

SFR1 3.1 1.9 1.6 

SFR2 8.0 2.4 1.8 

SFR3 10.5 4.8 2.7 



threshold and one on fluence thresho ld (im a ge mo de). The 
post-launch analysis can be found in iBandl (2006) . We use 
his results to define a detection probability at peak flux 
f is- 150 kcv, based on the observed GRB duration distribu- 
tion. The mean flux threshold is ~ 0.2 ph cm -2 s _1 . We also 
defined a sub-sample of "bright SWIFT bursts" by selecting 
bursts with peak flux P15-150 kev > 1 ph cm -2 s _1 . 

3.2 Observational constraints 

Depending on the assumption for the intrinsic peak energy 
distribution, our model has four (Amati-like relation) or five 
(log-normal E p distribution) free parameters. These param- 
eters can be constrained by the following observations: (1) 
the log N — log P diagram of B ATSE bursts Jk ommers et alJ 
l2000l:lstern et al]l2000l, 12002ft . (2) the observed peak energy 
distribution of long bright GRBs dPreece et all2000l) (3) the 
observed fraction of XRRs and XRFs in the sample of GRBs 
detected by HETE2. We adopt the same definitions as the 
HETE2 team f or these sub-classes, as well as their published 
observed ratio iSakamoto et alJl2005l) . 

Combining those three observations, we get 41 data 
points, so that our model has 37 (resp. 36) degrees of free- 
dom in the Amati-like relation scenario (resp. the log-normal 
E p distribution scenario). We find the best fit parameters 
and their dispersion by \ 2 minimization, as detailed below. 

3.3 Parameter determination 

We consider the six cases corresponding to our three as- 
sumed SFRs and our two possible spectral scenarios. In each 
case, we are able to find a good fit to the data. The search for 
the best fit is made by randomly choosing I/ m in (resp. Z/ m ax) 
in the interval 10 46 -10 52 erg s _1 (resp. 10 51 -10 56 erg s _1 ), 
8 in the range 0.5-3 and by adjusting the normalization k 
to minimize the \ 2 . More than 10 5 sets of parameters are 
tried for each of the six scenarios, and more than 10 5 GRBs 
are simulated in the Monte Carlo run for each of these sets. 
For each case, we define "best models" as models that fit the 
three observational constraints at the 1 sigma level, i.e. mod- 
els for which xLn < X 2 < Xmin + A X 2 , with A* 2 = 40.5 
(resp. 39.5) for 37 degrees of freedom (resp. 36). We can 
then compute the mean value and the dispersion of the best 
model parameters, as well as of any function of these pa- 
rameters, such as the predicted SWIFT GRB redshift dis- 
tribution. An example is given in Fig. [5] in the case SFR3 



Table 3. Best models: predicted fraction of all and 
bright SWIFT GRBs above 2 = 6 and z = 7. The re- 
ported values are averaged over all "best models" . The uncer- 
tainty in each case can be estimated from Fig. [3] 



SWIFT Bright SWIFT 

Rate %(z > 6) %(z > 7) %(z > 6) %(z > 7) 





Amati-like relation E p oc L° 43 




SFR1 


0.7% 0.4% 0.3% 


0.1% 


SFR2 


2.5% 1.6% 0.8% 


0.4% 


SFR3 


15% 12% 2.0% 


1.8% 



log-normal peak energy distribution 



SFR1 1.4% 0.7% 0.6% 0.3% 

SFR2 4.5% 2.9% 1.6% 1.0% 

SFR3 21% 17% 6.2% 4.5% 



+ Amati-like relation. Very similar plots are obtained in all 
the other considered cases. In the top-left panel of Fig. [5] 
we show the location of the best models in the parameter 
space I/min-imax-^, using the logiV — logP constraint only. 
It clearly appears that this constraint can fix the slope 8 of 
the luminosity function but leaves a degeneracy on L m i n and 
imax- Adding the two other constraints improve the param- 
eter determination (top-right, bottom- left panels in Fig.|5J|. 
Indeed, it can be seen that L max is better constrained using 
the peak energy distribution of bright BATSE bursts. The 
uncertainty on Z/ m i n is slightly reduced using the third con- 
straint. Finally the comparison between the top-left panel 
(log N — log P only) and the bottom-right panel (three con- 
straints together) shows the improvement in the parameter 
determination given by our method. 

In all considered cases, we always find a clear minimum 
for x 2 ■ As expected, we also find that for each physical quan- 
tity, the mean value for the best models defined as above is 
very close to the value for the minimum \ 2 model. The qual- 
ity of the best fit can be seen in Fig.|3] Finally the results are 
shown in Table 1, where we indicate for each case the mean 
value and the dispersion of the best model parameters. The 
obtained error bars confirm that the slope of the luminosity 
function remains better determined that the minimum and 
maximum luminosities. 



4 RESULTS AND DISCUSSION 
4.1 Redshift distribution 

For each scenario, we report in Table 2 the mean red- 
shift we obtain for all GRBs (intrinsic mean redshift), all 
GRBs detected by SWIFT and bright GRBs detected by 
SWIFT. The corresponding redshift distributions are plot- 
ted in Fig. as well as the ob served distribution for SWIFT 
bursts (IJakobsson et al]l2006l) . Each simulated redshift dis- 
tribution is computed averaging over the "best models" . The 
computed dispersion around this mean is the hatched region 
in Fig.H 

The observed redshift distribution of SWIFT bursts is 
plagued by selection effects, difficult to account for. For ex- 
ample, the deficiency of low redshift (2 ~ 1) events, with 
respect to pre-SWIFT bursts (see fig. 3 in IJakobsson et alJ 



6 F. Daigne, E.M. Rossi and R. Mochkovitch 




Figure 3. Best model (SFR3, Amati-like relation) . Left: the simulated log N — logP diagram of BATSE bursts is plotted as well 
as BATSE data iStern et al 1l200dl200^ : right: the simulated peak energy distribution of bright BATSE bursts is plotted as well as the 
observed distribution I Prcccc ct al. l200(t . 



2006) is not yet understood. For this reason, we do not 
attempt a formal fit of the data. However, we expect the 
observed distribution to be located between the simulated 
distributions of all SWIFT bursts (thin hatched region in 
fig-EJ an d of "bright" SWIFT bursts, with peak flux above 
1 ph cm -2 s 1 in the 15-150 keV band (thick hatched region). 
Figure 0] clearly shows that the expected redshift distri- 
bution strongly depends on the assumption made on the 
GRB comoving rate, so that SWIFT data can already put 
severe constraints on this rate. The case SFRi illustrates 
the effect of instrument thresholds : with a mean redshift 
{2} ~ 1.3 — 1.9, the simulated redshift distribution of bursts 
detected b y SWIFT is very far from the observed one 
((2) ~ 2.8; Jakobss on et al]l2006l) . despite an intrinsic dis- 
tribution of all GRBs that is much closer. Similarly SFR2 
gives a lower mean redshift ((z) ~ 1.5 — 2.4, see table 2) 
than the observed one. Conversely, SFR3 gives a much bet- 
ter agreement with data ((2) ~ 2.1 — 4.8). The observed 
SWIFT distribution lies just between the simulated redshift 
distribution of all SWIFT bursts and that of bright SWIFT 
bursts, as expected. This is qualitatively in agre ement with 
the analysis performed by Jak obsson et al] <l2006l) who found 
that SWIFT data were better reproduced using Model II of 
iNataraian et all (I2005T) . which corresponds to a SFR still 
rising up to 2 ~ 7, i.e. close to our SFR3. Our results seem 
rather independent on the details of the intrinsic GRB prop- 
erties, as the Amati-like relation and the log-normal E p give 
very similar distributions. 

We provide in addition in Table 3 the predicted fraction 
of all and bright SWIFT bursts above z — 6 and z = 7. We 
find that this fraction is less than 1 % for SFRi. Assuming 
100 SWIFT GRBs per year, we expect less than 5 GRBs 
per year above 2 = 6 for SFR2 (resp. less than 2 bright 



GRBs per year). In this case, most high redshift GRBs 
will be located around z = 6 — 7. In the SFR3 case, we 
expect 15-20 SWIFT GRBs per year above z = 6, including 
2-6 bright GRBs. In this case, a non negligible fraction of 
GRBs is predicted above z = 7. Clearly, the detection of 
GRB 050904 at z = 6.29 during the first year of SWIFT 
jKawai et alj feooe) is very unlikely with SFRi, marginally 
compatible with SFR2 and more easily explained with SFR3. 

The fact that SFR3 seems preferred by data is quite sur- 
prising, as this SFR is probably unrealistic: such a high rate 
of star formation at large redshift would result in an intense 
production of metals that would lead to much higher metal- 
lici ties in the structures and in the IGM than observed (see 
e.g. lDaigne et al]|2004lE005r) . Therefore, our results provide 
strong evidence that the properties of GRBs or/and GRB- 
progenitors are redshift dependent. 

One possibility is that the efficiency of GRB production by 
stars decreases with time. The ratio k of GRBs over type 
II supernovae would thus increase with redshift. If the true 
cosmic SFR in the Universe is more similar to SFRi (resp. 
SFR2), then our results show that this efficiency at z ~ 6 — 7 
is ~ 6 — 9 (resp. ~ 2 — 3) higher than at 2 ~ 2 (see Fig. 0. 
There is indeed increasing observational e vidence that the 
GRB rate does not directly trace the SFR llLe Floc'h et al] 
2006). Such an evolution could be related to many factors : 
metallicity, mass, rotation, or binarity of the progenitors. 
On the other hand, there could be an evolution of GRB in- 
trinsic properties with time (e.g. more luminous GRBs at 
large redshift). We tested the impact of a GRB luminosity 
evolution by considering a model where L m i n and I/ ma x scale 
as (I + 2)", the slope of the luminosity function being kept 
constant. We found that it is possible to reconcile SFRi with 
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Figure 4. Best models: redshift distribution. Each diagram shows for the six considered cases (1) the observed distribution of 
SWIFT bursts with measured redshifts (thick line); (2) the redshift distribution of all simulated bursts (thin dotted line); (3) the redshift 
distribution of simulated bursts that are detected by SWIFT (thin hatched region); (4) the redshift distribution of bright SWIFT bursts 
(thick hatched region). 



SWIFT data but this would imply a very strong evolution 
(i> > 2), probably unrealistic. However, it is unlikely that 
progenitor characteristics would evolve without any chang- 
ing in the GRB properties. Therefore, we suspect that both 
evolutionary effects are, in fact, present but that our results 
are mostly due to the former. 



4.2 Luminosity function and peak energy 
distribution 



The Amati-like relation scenario and the log-normal E p dis- 
tribution scenario lead to rather similar redshift distribution 
but differ strongly in some other aspects: (i) as can be seen 
in Fig. (right panel) , selection effects only can not create 
an apparent Amati-like correlation between the luminosity 
and the intrinsic peak energy when it is not originally 
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Figure 5. Amati-like relation in the best model (SFR3). GRB density contours are plotted in the lu minosity-intrinsic peak energy 
plane for SWIFT bursts. The slope is indicated in solid line, as well as the observed slope in dotted line JOhirlanda et af1l200Rf> . Left : 
simulation using an intrinsic Amati-like relation; right : simulation where the intrinsic peak energy follows a log-normal distribution. 
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Figure 6. Observed peak energy distribution in the best model (SFR.3). The peak energy distribution of all simulated GRBs, 
and of simulated GRBs that are detected by BATSE, HETE2 and SWIFT are plotted, insets showing with more details the regions of 
interest. Left : simulation using an Amati-like relation; right : simulation where the peak energy follows a log-normal distribution. 
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present. On the other hand, when an intrinsic correlation is 
assumed, detection thresholds do not modify the observed 
slope (left panel). This means that if the Amati relation 
is confirmed, it should have an intrinsic origin that should 
be explained by GRB models; (ii) Figure [(J shows that 
the intrinsic peak energy distribution is very different in 
the two scenarios. If an Amati-like relation is present, the 
intrinsic peak energy distribution is centered at about a 
few keV (4.5 keV for SFR3), since dim bursts are more 
numerous. Thus, the observed peak energy distribution 
in th e sub-sample of bright BATSE GRBs jPreece et alJ 
2000) is not representative of the whole GRB population, 
which is largely dominated by low E p events. On the other 
hand, if the peak energy distribution is not correlated with 
luminosity and has a log-normal distribution, then the 
observed distribution is close to the intrinsic distribution, 
centered at ~ 100 keV. Notice that in both cases, the mean 
observed peak energy by HETE2 and SWIFT is smaller 
than by BATSE. This effect has also consequences on the 
luminosity function. As can be seen in Fig. |7| the intrinsic 
Amati-like relation favors the detection of low luminosity 
bursts, as they have a smaller peak energy and therefore 
more photons in the 15-150 keV band. In this scenario, 
the lower energy threshold (15 keV vs 25 keV) of SWIFT 
compared to BATSE allows SWIFT to be more efficient 
in sampling the luminosity function. In the log-normal E p 
distribution scenario, the luminosity function of SWIFT 
and BATSE bursts are very similar. Notice that in both 
scenarios, the observed luminosity function is far from the 
intrinsic distribution, due to the small detection efficiency 
for the dimmest (most numerous) GRBs. Finally note that 
in the log-normal distribution case, the intrinsic mean 
peak energy is found to be larger than the constant value 
assumed by PM01. 

To check the dependance of our results on the assumed 
shape of the luminosity function, we also made a few simu- 
lations using broken power-laws : 

ia f L-°- s for L min < L < L„ 
y L~ s for L* < L < L max 

The low-luminosity slope has been fixed to avoid an ad- 
ditional free parameter. T he value —0.5 is a prediction of 
the internal shock model jDaigne et al.ll2006l) . It happens 
that the low-luminosity branch remains mostly undetected 
so that the best fit parameters are very close to those ob- 
tained with a single power-law (L* adjusting itself to the 
previous value of L m in). Therefore, current data do not al- 
low to distinguish between single and broken power-laws for 
the GRB luminosity function. 



5 CONCLUSION 

Using Monte Carlo simulations of the long GRB population 
under the assumptions that (i) the GRB rate follows the 
SFR; (ii) the GRB luminosity function is a power-law; (iii) 
the peak-energy is either determined by an Amati-like in- 
trinsic relation or by a log-normal distribution, we obtained 
the following results : 

(i) Luminosity function: the slope of the power-law is 



well constrained, S ~ 1.5 — 1.7, while the minimum and 
maximum luminosities are not so well determined : L m i n ~ 
0.8 - 3 x 10 50 erg s" 1 and L max ~3- 5 x 10 53 erg s _1 . The 
observed distribution is very biased as most low-luminosity 
GRBs are not detected by current instruments. 

(ii) Intrinsic peak energy distribution: if the ob- 
served Amati relation is confirmed, the relation needs to 
be intrinsic. A consequence of the relation is that the ob- 
served peak energy distribution is not representative of the 
whole GRB population, which has a lower mean value of a 
few keV. Therefore, there should be many undetected XRRs 
and XRFs. On the other hand, a log-normal peak energy 
distribution can also be in very good agreement with all the 
observational constraints, but selection effects alone can not 
produce an observed Amati relation. In this case, the mean 
peak energy of the whole GRB population is close to 100 keV 
and the observed distribution is much more representative. 

(iii) GRB com oving rate: in agreement with 
iPorciani fc Madaul (1200 lb . we find that one GRB pointing 
towards us is produced every 10 5 to 10 6 supernovae in the 
Universe. This rate should be corrected for beaming to get 
the true GRB rate in the Universe. This would require 
to assume a distribution of beaming angle or to derive it 
using the observed correlation with the isot ropic luminosity 
jFrail et all 1200 ll: iGhirlanda et al.l |2004) . The present 
redshift distribution of SWIFT bursts strongly favors 
SFR3, which still increases above z ~ 2. However, SFR3 
is probably unrealistic, as so many stars at high redshift 
would over-produce metals. This leads to conclude that 
some GRBs or /and GRB progenitor properties evolve with 
redshift. Our analysis suggests that the main evolutionary 
effect could be the redshift dependence of the efficiency 
of GRB production by stars, so that the GRB comoving 
rate still increases above z ~ 2, even if the SFR flattens or 
decreases. Thus, GRBs would not directly trace the SFR, as 
also suggested by rece nt studies on the prop erties of GRB 
host galaxies (see e.g. iLe Floc'h et al.l I2006T) . To reconcile 
our results with a more plausible SFR such as SFRi or even 
SFR2, this efficiency should be about 6 — 7 times larger at 
z = 7 than at z = 2. 

(iv) Detection rate of high redshift GRBs: based on 
a yearly rate of ~ 100 SWIFT GRBs, we predict ~ 2 - 6 
(resp. ~ 1 — 5) bright SWIFT bursts per year above z — 6 
(resp. z = 7). 
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